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Abstract 



A measurement of the exclusive two-photon production of muon pairs in proton- 
proton collisions at ^/s = 7 TeV, pp — > p^+/^^p, is reported using data correspond- 
ing to an integrated luminosity of 40 pb^^. For muon pairs with invariant mass 
greater than 11.5 GeV, transverse momentum Pt(^) > 4 GeV and pseudorapidity 
< 2.1, a fit to the dimuon pj{ji^ii^) distribution results in a measured cross 
section of o-{p -> p^+|^-p) = 3.38+^;| (stat.) ± 0.16 (syst.) ± 0.14 (lumi.) pb, con- 
sistent with the theoretical prediction evaluated with the event generator Lpair. The 
ratio to the predicted cross section is 0.83+g;{3 (stat.) ± 0.04 (syst.) ± 0.03 (lumi.). The 
characteristic distributions of the muon pairs produced via 77 fusion, such as the 
muon acoplanarity, the muon pair invariant mass and transverse momentum agree 
with those from the theory. 
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1 Introduction 

The exclusive two-photon production of lepton pairs may be reUably calculated within the 
framework of quantum electrodynamics (QED) [l] (Fig. [T]), within uncertainties of less than 
1% associated with the proton form factor f2]. Indeed, detailed theoretical studies have shown 
that corrections due to hadronic interactions between the elastically scattered protons are well 
below 1% and can be safely neglected [3J. The unique features of this process, like the extremely 
small pair transverse momentum and acoplanarity (defined as 1 — \ A(p{ji^ji^) / ti\), stem from 
the very small virtualities of the exchanged photons. 

At the Tevatron, the exclusive two-photon production of electron fl'.'S] and muon [5','6l pairs in 
pp collisions has been measured with the CDF detector. Observations have been made of QED 
signals, leading to measurements of exclusive charmonium photoproduction [6J and searches 
for anomalous high-mass exclusive dilepton production [5|. However, all such measurements 
have very limited numbers of selected events because the data samples were restricted to sin- 
gle interaction bunch crossings. The higher energies and increased luminosity available at the 
Large Hadron Collider (LHC) will allow significant improvements in these measurements, if 
this limitation can be avoided. As a result of the small theoretical uncertainties and character- 
istic kinematic distributions in 77 — >• , this process has been proposed as a candidate for 
a complementary absolute calibration of the luminosity of pp collisions 

Unless both outgoing protons are detected, the semi-exclusive two-photon production, involv- 
ing single or double proton dissociation (Fig. 1, middle and right panels), becomes an irre- 
ducible background that has to be subtracted. The proton-dissociation process is less well 
determined theoretically, and in particular requires significant corrections due to proton rescat- 
tering. This effect occurs when there are strong-interaction exchanges between the protons, in 
addition to the two-photon interaction. These extra contributions may alter the kinematic dis- 
tributions of the final-state muons, and may also produce additional low-momentum hadrons. 
As a result, the proton-dissociation process has significantly different kinematic distributions 
compared to the pure exclusive case, allowing an effective separation of the signal from this 
background. 




Figure 1: Schematic diagrams for the exclusive and semi-exclusive two-photon production of muon 
pairs in pp collisions for the elastic (left), single dissociative (center), and double dissociative (right) 
cases. The three lines in the final state of the center and right plots indicate dissociation of the proton 
into a low-mass system N. 

In this paper, we report a measurement of dimuon exclusive production in pp collisions at 
= 7TeV for the invariant mass of the pair above 11.5 GeV, with each muon having trans- 
verse momentum Pt{}i) > 4GeV and pseudorapidity |?/(/^)| < 2.1 (where rj is defined as 
— ln(tan(0/2))). This measurement is based on data collected by the Compact Muon Solenoid 
(CMS) experiment during the 2010 LHC run, including beam collisions with multiple interac- 
tions in the same bunch crossing (event pileup), and corresponding to an integrated luminosity 
of 40 pb^^ with a relative uncertainty of 4% [7J. 
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3 Simulated Samples 



The paper is organized as follows. In Section 2, a brief description of the CMS detector is 
provided. Section 3 describes the data and samples of simulated events used in the analysis. 
Section 4 documents the criteria used to select events, and Section 5 the method used to extract 
the signal yield from the data. The systematic uncertainties and cross-checks performed are 
discussed in Section 6, while Section 7 contains plots comparing the selected events in data and 
simulation. Finally, the results of the measurement are given in Section 8 and summarized in 
Section 9. 

2 The CMS detector 

A detailed description of the CMS experiment can be found elsewhere fS]. The central fea- 
ture of the CMS apparatus is a superconducting solenoid, of 6 m internal diameter. Within 
the field volume are the silicon pixel and strip tracker, the crystal electromagnetic calorime- 
ter, and the brass/ scintillator hadronic calorimeter. Muons are measured in gaseous detectors 
embedded in the iron return yoke. Besides the barrel and endcap detectors, CMS has exten- 
sive forward calorimetry. CMS uses a right-handed coordinate system, with the origin at the 
nominal collision point, the x axis pointing to the center of the LHC ring, the y axis pointing 
up (perpendicular to the plane of the LHC ring), and the z axis along the anticlockwise-beam 
direction. The azimuthal angle (p is measured in the x-y plane. Muons are measured in the 
window \rj\ < 2.4, with detection planes made using three systems: drift tubes, cathode strip 
chambers, and resistive plate chambers. Thanks to the strong magnetic field, 3.8 T, and to the 
high granularity of the silicon tracker (three layers consisting of 66 million 100 x 150 fim^ pix- 
els followed by ten microstrip layers, with strips of pitch between 80 and 180 fim), the px of 
the muons matched to silicon tracks is measured with a resolution better than ~ 1.5%, for pj 
less than 100 GeV. The first level of the CMS trigger system, composed of custom hardware 
processors, uses information from the calorimeters and muon detectors to select (in less than 
1 }is) the most interesting events. The High Level Trigger processor farm further decreases the 
event rate from 50-100 kHz to a few hundred Hz, before data storage. 

3 Simulated Samples 

The Lpair 4.0 event generator 19^. TOl is used to produce simulated samples of two-photon 
production of muon pairs. The generator uses full leading-order QED matrix elements, and 
the cross sections for the exclusive events depend on the proton electromagnetic form-factors 
to account for the distribution of charge within the proton. For proton dissociation, the cross 
sections depend on the proton structure function. In order to simulate the fragmentation of 
the dissociated proton into a low-mass system N, the LUND model shower routine [11] imple- 
mented in the JetSet software IIT2II is used with two different structure functions. For masses of 
the dissociating system m^f < 2 GeV and photon virtualities < 5 GeV^, the Brasse "cluster" 
fragmentation is chosen |13|, while for the other cases the Suri-Yenni "string" fragmentation 
is applied [14J. In the first case, the low-mass system N mostly decays to a A+ or A++ reso- 
nance, which results in low-multiplicity states. In the second case, the high-mass system usu- 
ally decays to a variety of resonances (A, p, CI, t], K), which produce a large number of forward 
protons, pions, neutrons and photons. No corrections are applied to account for rescattering 
effects. In general, these effects are expected to increase with the transverse momentum of the 
muon pair, modifying the slope of the pj(/i+f/^) distribution [3J. 

The inclusive Drell-Yan (DY) and quantum chromodynamic (QCD) dimuon backgrounds are 
simulated with pythia v. 6.422 ||T5|. using the Z2 underlying event tune ||T6|. All these samples 
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are then passed through the full GEANT4 detector simulation l|T7ll in order to determine the 
signal and background efficiencies after all selection criteria are applied. 

4 Event selection 

The analysis uses a sample of pp collisions at a/s = 7 TeV, collected during 2010 at the LHC 
and corresponding to an integrated luminosity of 40 pb^^. The sample includes 36 pb^^ of data 
passing the standard CMS quality criteria for all detector subsystems, and 4 pb^^ in which the 
quality criteria are satisfied for the tracking and muon systems used in the analysis. From the 
sample of triggered events, the presence of two reconstructed muons is required. Then the 
exclusivity selection is performed to keep only events with a vertex having no tracks other 
than those from the two muons. Finally, the signal muons are required to satisfy identification 
criteria, and kinematic constraints are imposed using their four-momentum. All selection steps 
are described in the following sections. 

4.1 Trigger and muon reconstruction 

Events are selected online by triggers requiring the presence of two muons with a minimum pj 
of 3 GeV. No requirement on the charge of the muons is applied at the trigger level. Muons are 
reconstructed offline by combining information from the muon chambers with that on charged- 
particle tracks reconstructed in the silicon tracker |ll8| , and events with a pair of oppositely 
charged muons are selected. 

4.2 Vertex and tracl< exclusivity selection 

With single interactions, the exclusive signal is characterized by the presence of two muons, no 
additional tracks, and no activity above the noise threshold in the calorimeters. The presence 
of additional interactions in the same bunch crossing will spoil this signature by producing 
additional tracks and energy deposits in the calorimeters. In the 2010 data, less than 20% of the 
total luminosity was estimated to have been collected from bunch crossings where only a single 
interaction look place, leading to a significant decrease in signal efficiency if the conditions of 
no extra tracks or calorimeter energy are required. 

The selection of exclusive events is therefore applied using the pixel and silicon tracker only, 
since the primary vertex reconstruction |T9l I20II allows discrimination between different inter- 
actions within the same bunch crossing. The selection requires a valid vertex, reconstructed 
using an adaptive vertex fit to charged-particle tracks clustered in z ll20ll2T| . with exactly two 
muons and no other associated tracks, and vertex fit probability greater than 0.1%. The dimuon 
vertex is further required to have coordinates consistent with a collision in CMS, with a longi- 
tudinal displacement of less than 24 cm and a transverse displacement of less than 0.1 cm. 

In order to reduce the background from inclusive DY and QCD dimuon production, which typ- 
ically have many tracks originating from the same vertex as a prompt muon pair, the dimuon 
vertex is required to be separated in three dimensions by more than 2 mm from any additional 
tracks in the event. This value is selected to optimize the signal efficiency and background re- 
jection found in events triggered only by the presence of colliding bunches ("zero-bias" events), 
and in DY Monte Carlo simulation. For the zero-bias data, this is accomplished by introducing 
an artificial additional dimuon vertex into each event as a proxy for an exclusive dimuon inter- 
action. Thus, in this study, beam crossings with no real vertex present are counted as "single 
vertex" events, and crossings with one real vertex are counted as having an additional pileup 
event. 
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4 Event selection 



The effects of the track veto on the signal efficiency and on the efficiency for misidentifying 
background as signal are studied as a function of the distance to the closest track for the zero- 
bias sample and DY background (Figs. |2] and |3|. With no extra vertices in the zero-bias events, 
the efficiency approaches 100% as expected for events with no pileup. With the addition of 
overlap events, the efficiency decreases, falling to ~ 60% with 8 extra vertices reconstructed. 
In the full data sample the average number of extra vertices is 2.1, with less than 10% of events 
having 4 or more extra vertices. The efficiency for misidentifying the DY background as signal 
increases sharply for distances less than 1 mm, consistent with the resolution of the single-track 
impact parameter in the longitudinal direction pO | . 
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Figure 2: Efficiency of the zero extra tracks selection vs. distance to closest track computed with the 
artificial vertex method in zero-bias data. The points correspond to events with 0, 1, 2, and 8 real vertices 
in the event. Events to the right of the vertical dashed line are selected. The vertical error bars are 
negligible. 



4.3 Muon identification 

Each muon of the pair is required to pass a "tight" muon selection |'22'|. This selection consists of 
requesting that the reconstructed muon have at least one hit in the pixel detector, at least 10 hits 
in the silicon strip tracker, and segments reconstructed in at least two muon detection planes. 
In addition, a global fit to the combined information from the tracker and muon systems must 
include at least one muon chamber hit, and have a per degree of freedom of less than 10. 

4.4 Kinematic selection 

In order to minimize the systematic uncertainties related to the knowledge of the low-pj and 
large-// muon efficiencies, only muons with -pj > A GeV and \r]\ < 2.1 are selected. The pj and 
\r]\ requirements retain muon pairs from exclusive photoproduction of upsilon mesons, 7p — > 
Yp — > n^]i^p. This process occurs when a photon emitted from one proton fluctuates into a qq 
pair, which interacts with the second proton via a color-singlet exchange. This contribution is 
removed by requiring that the muons have an invariant mass m{}i^}i^ ) > 11.5 GeV. 
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Figure 3: Efficiency of the zero extra tracks selection vs. distance to closest track computed for DY 
events in simulation. Events to the right of the vertical dashed line are selected. The vertical error bars 
are negligible. 

In order to suppress further the proton dissociation background, the muon pair is required 
to be back-to-back in azimuthal angle (1 — \A(p{}i^}i^) / 7t\ < 0.1) and balanced in the scalar 
difference in the pj of the two muons (|Apx(?^^/^ )| < 1-0 GeV). A possible contamination 
could arise from cosmic-ray muons, which would produce a signature similar to the exclusive 
77 — > ji^ji'^ signal. The three-dimensional opening angle of the pair, defined as the arccosine 
of the normalized scalar product of the muon momentum vectors, is therefore required to be 
smaller than 0.95 n, to reduce any contribution from cosmic-ray muons. 

The effect of each step of the selection on the data and simulated signal and background sam- 
ples is shown in Table [T] After all selection criteria are applied, 148 events remain, where from 
simulation, approximately half are expected to originate from elastic production. The number 
of events selected in data is below the expectation from simulation, with an observed yield that 
is roughly 80% of the sum of simulated signal and background processes. The deficit could be 
caused by a lower-than-expected signal yield, or by a smaller proton-dissociation contribution 
than expected from simulation. 

5 Signal extraction 
5.1 Efficiency corrections 

A correction is applied to account for the presence of extra proton-proton interactions in the 
same bunch crossing as a signal event. These pileup interactions will result in an inefficiency 
if they produce a track with a position within the nominal 2 mm veto distance around the 
dimuon vertex. This effect is studied in zero-bias data using the method described in Sec- 
tion 4.2. The nominal 2 mm veto is then applied around the dimuon vertex, and the event is 
accepted if no tracks fall within the veto distance. The efficiency is measured as a function of 
the instantaneous luminosity per colliding bunch. The average efficiency is calculated based on 
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5 Signal extraction 



Table 1: Number of events selected in data and number of signal and background events ex- 
pected from simulation at each selection step for an integrated luminosity of 40pb^^. The last 
column is the number of events expected from the sum of the signal, DY, and proton dissocia- 
tion backgrounds in the simulation. The relative statistical uncertainty on the sum of simulated 
signal and background samples in each row is < 0.5%. The contribution from exclusive res- 
onance production of Y or Xb mesons is not simulated, and thus contributes only to the data 
column before requiring m{}i^}i^) > 11.5 GeV. For entries in the line "Muon ID" and below, 
the simulation is corrected for effects related to event pileup, muon identification, trigger, and 
tracking efficiencies, as described in the text. 



Selection 


Data 


Signal 


Single-pdiss. 


Double-pdiss. 


DY 


Total 


Vertex and track-exclusivity 


921 


247 


437 


197 


56 


937 


Muon ID 


724 


193 


336 


160 


53 


741 


Pt > 4GeV, \ f]\ < 2.1 


438 


132 


241 


106 


20 


499 


m{}i+}i-) > 11.5 GeV 


270 


95 


187 


86 


13 


380 


3D angle < 0.95/1 


257 


87 


178 


83 


12 


361 


1 - \A(p/n\ < 0.1 


203 


87 


126 


41 


8 


263 


lApxl < 1.0 GeV 


148 


86 


79 


16 


3 


184 



the instantaneous luminosities to be 92.29% for the full 2010 data set, with negligible statistical 
uncertainty. 

The trigger, tracking, and offline muon selection efficiencies are each obtained from the tag- 
and-probe ||22ll23l method by using samples of inclusive ]/xp /^^H F^l^^ events 

from data and Monte Carlo simulation. These control samples are triggered on one muon such 
that the other muon is unbiased with respect to the efficiency to be measured. For pj < 20 GeV 
muons from ]/xp decays are used, while above 20 GeV muons from Z decays are used. The 
trigger and offline muon selection efficiencies are measured using ]/ip events by requiring a 
muon tag that, when combined with a track reconstructed using only the silicon detectors, is 
consistent with a j/ip. These efficiencies are measured in bins in pj and t], separately for the 
two muon charges. The tracking efficiency is measured similarly using ]/xp events by requiring 
a muon tag that, when combined with a muon reconstructed using only the muon systems, is 
consistent with a ]/xp. The tracking efficiency is then measured on the unbiased muon probe. In 
contrast to the trigger and muon identification efficiencies, the tracking efficiency is measured 
in data and Monte Carlo simulation averaged over \ri\ < 2.1 and pj > A GeV, and is taken to be 
uncorrelated between the two tracks. The resulting ratio of efficiencies in data and simulation 
for the pair (99.18 ± 0.14)% is applied as a correction to the efficiency. 

The effect of the vertexing efficiency is studied both in inclusive dimuon data and signal simu- 
lation, by performing an independent selection of all muon pairs with a longitudinal separation 
of less than 0.5 mm. A Kalman filter [24J algorithm is then applied to estimate the best position 
of the dimuon vertex, without using information from any tracks other than the two muons. 
Among events with a valid dimuon vertex and for which no additional tracks exist within 
2 mm in z, the efficiency for the default adaptive vertex fitter to reconstruct a primary vertex 
with only two muons attached and matching with the Kalman vertex is computed. The ratio 
of the vertexing efficiency in data to that in simulation is 99.97%, and therefore no correction is 
applied. 



5.2 Maximum likelihood fit 
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Figure 4: Distribution of pj{}i^ii~) for the selected sample. Data are shown as points with sta- 
tistical error bars. The histograms represent the simulated signal (yellow), single (light green) 
and double (dark green) proton dissociative backgrounds, and DY (red). The yields are deter- 
mined from a fit using the distributions from simulation. 

5.2 Maximum likelihood fit 

The elastic pp — )■ pfi^fi^p contribution is extracted by performing a binned maximum-likelihood 
fit to the measured px(/^^F ) distribution. Shapes from Monte Carlo simulation are used for 
the signal, single-proton dissociation, double-proton dissociation, and DY contributions, with 
all corrections described in Section 4.4 applied. 

Three parameters are determined from the fit: the elastic signal yield relative to the Lpair pre- 
diction for an integrated luminosity of 40 pb^^ {Rei-ei), the single-proton dissociation yield 
relative to the Lpair single-proton dissociation prediction for 40 pb^^ (Rrf/ss-E/)/ arid an expo- 
nential modification factor for the shape of the pj distribution, characterized by the parame- 
ter a. The modification parameter is included to account for possible rescattering effects not 
included in the simulation, as described in Section 3. Given the small number of events ex- 
pected in 40 pb^^, the double-proton dissociation and DY contributions cannot be treated as 
free parameters and are fixed from simulation to their predicted values. The contribution from 
exclusive 77 — ?► t+t^ production is estimated to be 0.1 events from the simulation, and is 
neglected. 

The pxlfi^H ) distribution in data is shown overlaid with the result of the fit to the shapes 
from Monte Carlo simulation in Fig. |4j The result from the best fit to the data is: 



data-theory signal ratio: 
single-proton dissociation yield ratio: 
modification parameter: 




(1) 



with asymmetric statistical uncertainties computed using MINOS [25 J. The corresponding value 
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5 Signal extraction 



of the signal cross section is S.SS^qI^ (stat.) pb. The resulting la and 2a contours projected onto 
each pair of fit variables are displayed in Fig. |5j For any values of the proton dissociation ratio 
and slope within the la contour, the extreme values of the data-theory signal ratio are 0.64 
and 1.03. The upper value of 1.03 for the signal ratio would correspond to the single-proton 
dissociation component having a ratio to the prediction of approximately 0.65. The best fit 
does not require a significant modification parameter. However, the statistical uncertainty on 
this parameter is chosen to play the role of a non-negligible systematic uncertainty, to take 
account of the neglect of the rescattering effects in the simulation. 
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Figure 5: One and two standard-deviation contours in the plane of fitted parameters for the 
proton-dissociation yield ratio vs. modification parameter a (left), the data-theory signal ratio 
vs. modification parameter a (center), and the data-theory signal ratio vs. proton-dissociation 
yield ratio (right). The contours represent 39.3% and 86.5% confidence regions, where the cross 
indicates the best-fit point. 

As a cross-check, a fit to the 1 — \ A(p{}i^}i^) / k\ distribution is performed, with the signal and 
single-proton dissociation yields as free parameters, and the shape of the single-proton dissoci- 
ation component fixed from the simulation. The resulting value of the data-theory signal ratio 
is 0.81 ^qIj, consistent with the nominal fit result. 

The central values of the signal and single-proton dissociation yields from the fit are both below 
the mean number expected for 40 pb^^, consistent with the deficit shown in Table 1. This is in- 
vestigated by repeating the fit, first with the Apx(/^^?^ ) < 1-0 GeV requirement removed, and 
then with both the Apj{pi^}i^) < l.OGeV and 1 — |A^(^+f/^)/7r| < 0.1 selections removed. 
From simulation this is expected to have negligible effect on the signal efficiency, while en- 
hancing the background. The double-proton dissociation and DY contributions in particular 
are expected to be small with the nominal selection, but their sum becomes comparable in size 
to the signal with the ApT(f<^H ) 1 ~ l^'PiH^'H ) / ^1 requirements removed. 

The fits to the data with these looser selection requirements are shown in Fig. |6| and the re- 
sulting best-fit yields for the signal and single-proton dissociation are shown in Table |2| the 
single-proton dissociation yield is observed to be significantly lower relative to the prediction 
with the looser selections. In all variations, the normalizations of the double-proton dissocia- 
tion and DY yields are fixed, although the double-dissociation contribution is expected to be 
significant at large pj{}i^ji^) with the looser selection. With additional data, a more precise 
comparison of the single and double dissociation yields to the theoretical expectation may be 
made. In spite of the lower single-proton dissociation yield, the data-theory ratio for the signal 
is stable in all three variations. 
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Figure 6: Result of fit to the pj{ii^}i^) distribution with requirements on |Apt(^+^^)| (left) 
and on both |Apx(f< ) I and 1 — |A^(^+f/^)/ n\ (right) removed. The points with error bars 
represent the data. The histograms are the result of fitting the simulated distributions to the 
data. 

Table 2: Best-fit values of Rei-ei arid Rdiss-Ei for the nominal selection, and with the require- 
ments on jApxlM^f )l arid 1 — |A(^(f/+/i^)/7r| removed. 



Selection 


R-El-El 


Rdiss-El 


All selection criteria applied 




73+0-16 


No Apx requirement 




63+°-" 


Both |ApT and 1 — A^/ n\ requirements removed 




45+0.08 



6 Control plots 

The dimuon invariant mass and acoplanarity distributions for events passing all selection crite- 
ria listed in Table[l]are shown in Fig.[^ with the simulation predictions normalized to the best- 
fit signal and background yields. The event with the largest invariant mass has m(f/+f/^) = 
76 GeV. No events consistent with Z — ?► are observed. This is expected for exclusive 

production, since the 77 — > Z process is forbidden at tree-level. 

In Fig. [sj the |Apx(F^/^ )l and rj distributions are plotted. In Figs. |9- 10 the data and sim- 
ulation are similarly compared for the pj and t] of single muons passing all other selection 
requirements. Agreement between the data and simulation is observed in the distributions of 
all dimuon and single-muon quantities. 



7 Systematic uncertainties and cross-checks 

Systematic uncertainties related to the pileup efficiency correction, muon trigger and recon- 
struction efficiency corrections, momentum scale, LHC crossing angle, and description of the 
backgrounds in the fit are considered. The systematic uncertainties related to the muon iden- 
tification, trigger, and tracking efficiencies are determined from the statistical uncertainties of 
the ]/xp and Z control samples used to derive the corrections. The remaining systematic uncer- 
tainties are evaluated by varying each contribution as described in the following sections, and 
repeating the fit with the same three free parameters Rei-eu Russ-eu and the shape correction 
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7 Systematic uncertainties and cross-checl<s 
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Figure 7: Muon pair invariant mass spectrum (left) and acoplanarity (right), with all selection 
criteria applied and the simulation normalized to the best-fit value. Data are shown as points 
with statistical error bars, while the histograms represent the simulated signal (yellow), single 
(light green) and double (dark green) proton dissociative backgrounds, and DY (red). 
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Figure 8: Muon pair transverse momentum difference (left) and pair pseudorapidity (right), 
with all selection criteria applied and the simulation normalized to the best-fit value. Data are 
shown as points with statistical error bars, while the histograms represent the simulated signal 
(yellow), single (light green) and double (dark green) proton dissociative backgrounds, and DY 
(red). 

a. The relative difference of the data-theory signal ratio between the modified and the nominal 
fit result is taken as a systematic uncertainty. 

7.1 Pileup correction systematic uncertainties 

Charged tracks from pileup interactions more than 2.0 mm from the dimuon vertex may induce 
a signal inefficiency, if they are misreconstructed to originate from within the 2.0 mm veto win- 
dow. The ?/-dependent single-track impact parameter resolution in CMS has been measured to 
be less than 0.2 mm in the transverse direction, and less than 1.0 mm in the longitudinal direc- 
tion II20II . The track-veto efficiency is studied in zero-bias data by varying the nominal 2.0 mm 
veto distance from 1.0 to 3.0 mm. The maximum relative variation is found to be 3.6%, when 
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Figure 9: Single-muon pseudorapidity distribution with all other selections applied for f/+ (left) 
and ji^ (right) and the simulation normalized to the best-fit value. Data are shown as points 
with statistical error bars, while the histograms represent the simulated signal (yellow), single 
(light green) and double (dark green) proton dissociative backgrounds, and DY (red). 




P,(^") [GeV] p^(^-) [GeV] 

Figure 10: Single-muon transverse momentum with all other selections applied for ^+ (left) 
and f/^ (right) and the simulation normalized to the best-fit value. Data are shown as points 
with statistical error bars, while the histograms represent the simulated signal (yellow), single 
(light green) and double (dark green) proton dissociative backgrounds, and DY (red). 

enlarging the veto size to 3.0 mm. In addition, the veto is modified to use high-quality tracks 
having at least seven consecutive layers hit in the tracker, in place of the default veto based on 
all charged tracks. This is found to result in a 2.5% variation in the signal yield, which is taken 
as a systematic uncertainty. 

As a further check, the same variations are applied to the selected sample of dimuon events, 
removing the Y mass cut m < 11.5 GeV to increase the statistics with photo-produced exclusive 
upsilon events. The change in the number of events selected in the dimuon sample is found to 
be consistent with the expectation from the zero-bias sample. 
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7 Systematic uncertainties and cross-checl<s 



7.2 Muon efficiencies and momentum scale 

The statistical uncertainty on the muon efficiency correction is evaluated by performing a 
fast Monte Carlo study in which each single-muon correction evaluated from the tag-and- 
probe study is varied independently using a Gaussian distribution having a width equal to 
the measured uncertainty. The r.m.s. of the distribution of the resulting variations in the over- 
all dimuon efficiency correction is taken as the systematic uncertainty. From 1000 pseudo- 
experiments, this results in an uncertainty of 0.8%. In addition, we study the effect of cor- 
relations in the dimuon efficiency. The tag-and-probe study is only sensitive to single-muon 
efficiencies. Since we take the dimuon efficiency as the product of the single-muon efficiencies, 
the effect of correlations in the efficiency are not modeled. To evaluate the size of this effect, the 
efficiency corrections are computed after removing events in the ]/xp control sample in which 
the two muons bend towards each other in the r-(p plane, potentially becoming very close or 
overlapping. Such events may introduce larger correlations in the efficiency of the dimuon pair 
than would be present in the well separated signal muons. Repeating the signal extraction with 
this change results in a relative difference of 0.7% from the nominal efficiency, which is taken 
as a systematic uncertainty. 

Using studies of the muon momentum scale derived from Z }^^}^~ f^S^I, the muon pj is 
shifted by the observed px^dependent bias, and the nominal fit is performed again. The result- 
ing relative change in the signal yield is 0.1%, which is taken as a systematic uncertainty. As a 
cross-check using a sample kinematically closer to the signal, we apply all the selections except 
for the veto on the Y mass region, and perform a fit to the Y(1S) resonance. The resulting mass 
is consistent with the PDG value [26], within an uncertainty of 20 MeV. Applying a 20 MeV 
shift to the mass and px scales of the data and performing the fit again results in no change 
from the nominal efficiency. 

7.3 Vertexing and tracl^ing efficiencies 

Since the study described in Section 4.4 shows no significant difference in the vertexing ef- 
ficiency between data and simulation, the 0.1% statistical uncertainty of the measurement in 
data is taken as a systematic uncertainty. For the tracking efficiency, the difference between 
data and simulation is applied as a single correction without binning in px or //. The statistical 
uncertainty of 0.1% on the correction for the dimuon is taken as a systematic uncertainty. 

7.4 Crossing angle 

The non-zero crossing angle of the LHC beams leads to a boost of the dimuon system in the 
X direction. Consequently, the px of the pair is over-estimated by a few MeV, especially for 
high-mass dimuon events. This effect is estimated by applying a correction for the Lorentz 
boost, using a half-angle of 100 ^rad in the x-z plane. This results in a 1.0% variation from the 
nominal fit value, and is taken as an additional systematic uncertainty. 

7.5 Fit stability 

Checks of the fit stability are performed by testing different bin widths and fit ranges. Starting 
from the nominal number of 20 bins in the range 0-3 GeV, variations in the bin width from 
0.1 to 0.2 GeV and fit range [0,2] to [0,4] GeV show deviations by at most 3.3% with respect 
to the nominal yield. The fit bias is studied by performing a series of Monte Carlo pseudo- 
experiments for different input values of the signal and proton-dissociation yields, using events 
drawn from the fully simulated samples. The means of the pull distributions are found to be 
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consistent with zero. Since the pseudo-experiments with the nominal binning and fit range 
show no significant bias, no additional systematics are assigned in this case. 

7.6 Backgrounds 

The yields of the double-proton dissociation and DY contributions are fixed in the nominal 
fit. To estimate the systematic uncertainty from this constraint, the fit is repeated with each 
of these varied independently by a factor of 2. The resulting changes in the fitted signal yield 
are 0.9% and 0.4%, respectively, where because of the similar shapes of the single and double 
proton dissociation components, this variation is partly absorbed into the fitted single-proton 
dissociation yield. As a cross-check of this procedure, the | Apx(f'^F ) I ^^id 1 — \ A(p{]i^}i^ ) /7r| 
requirements are inverted to select samples of events expected to be dominated by double- 
proton dissociation and DY backgrounds. The agreement between data and simulation in these 
regions is found to be within the factor of 2 used as a systematic variation. 

The possibility of a large contamination from cosmic-ray muons, which may fake a signal since 
they will not be correlated with other tracks in the event, is studied by comparing the vertex po- 
sition and three-dimensional opening angle in data and simulations of collision backgrounds. 
A total of three events fail the vertex position selection in data, after all other selection crite- 
ria are applied. All three also fail the opening angle selection, which is consistent with the 
expected signature from cosmic muons. We conclude that the opening angle requirement ef- 
fectively rejects cosmic muons, and do not assign a systematic uncertainty for this possible 
contamination. 

A similar check for contamination from beam-halo muons is performed by applying the nomi- 
nal analysis selection to non-collision events triggered by the presence of a single beam. Within 
the limited statistics, zero events pass all the analysis selections, and therefore no additional 
systematic uncertainty is assigned in this case. 

7.7 Summary of systematic uncertainties 

The individual variations in the definition of the track-veto are taken as correlated uncertain- 
ties, with the largest variation taken as a contribution to the systematic uncertainty. The largest 
variation related to the track quality, obtained when requiring high-purity tracks with > 10 hits 
instead of the nominal value of > 3 hits, is also taken as a contribution. The larger variation 
resulting from increasing or decreasing the double-proton dissociation background normaliza- 
tion by a factor of 2, and the larger variation resulting from increasing or decreasing the DY 
background normalization by a factor of 2, are each taken as contributions to the systematic 
uncertainty. The variation in the crossing angle, muon identification and trigger efficiencies, 
tracking efficiency, bias due to correlations in the ]/tp control sample, and vertexing efficiency 
are treated as uncorrelated uncertainties. Summing quadratically all uncorrelated contribu- 
tions gives an overall relative systematic uncertainty of 4.8% on the signal yield (Table |3|. 

8 Results 

For muon pairs with invariant mass greater than 11.5 GeV, single-muon transverse momentum 
Pt{}i) > 4GeV, and single-muon pseudorapidity in the range |?/(/^) | < 2.1, 148 events pass all 
selections. Approximately half of these are ascribed to fully exclusive (elastic) production. The 
number of events expected from Monte Carlo simulation of signal, proton dissociation, and DY 
backgrounds for an integrated luminosity of 40 pb^^ is 184. 
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9 Summary 



Table 3: Relative systematic uncertainties. 





A/ariafion from noinin;il vipIH 


-LlCl-V_-lN, V CLW V,lltdJ.Cl 


3 6% 


Track nualitv 


2.5% 


DY background 


0.4% 


Double-proton dissociation background 


0.9% 


Crossing angle 


1.0% 


Tracking efficiency 


0.1% 


Vertexing efficiency 


0.1% 


Momentum scale 


0.1% 


Efficiency correlations in ]/xp control sample 


0.7% 


Muon and trigger efficiency statistical uncertainty 


0.8% 


Total 


4.8% 



The resulting visible cross section from a fit to the pj ) distribution is cr(pp — )■ Y'V-^V P) = 
3.381q 55 (stat.) ± 0.16 (syst.) ± 0.14 (lumi.) pb, and the corresponding data-theory signal ratio 
is 0.83^Qj3 (stat.) ± 0.04 (syst.) ± 0.03 (lumi.), where the statistical imcertainties are strongly 
correlated with the single-proton dissociation background. 

9 Summary 

A measurement is reported of the exclusive two-photon production of muon pairs, pp — > 
p^+/^ p, in a 40 pb^ sample of proton-proton collisions collected at ^Js = 7TeV during 2010 
at the LHC. The measured cross section 



cr(pp p^+^"p) = 3.38lg|^ (stat.) ± 0.16 (syst.) ± 0.14 (lumi.) pb, 

is consistent with the predicted value, and the characteristic distributions of the muon pairs 
produced via 77 fusion, such as the pair acoplanarity and transverse momentum, are well 
described by the full simulation using the matrix-element event generator LPAIR. The detec- 
tion efficiencies are determined from control samples in data, including corrections for the 
significant event pileup. The signal yield is correlated with the dominant background from 
two-photon production with proton dissociation, for which the current estimate from a fit to 
the Pt(m^?^ ) distribution can be improved with additional data. The efficiency for the exclu- 
sivity selection is above 90% in the full data sample collected by CMS during the 2010 LHC 
run. With increasing instantaneous luminosity this efficiency will decrease, but without possi- 
ble improvements to the selection remains above 60% with up to 8 additional pileup vertices. 
Since the process may be calculated reliably in the framework of QED, within uncertainties 
associated with the proton form factor, this represents a first step towards a complementary 
luminosity measurement, and a reference for other exclusive production measurements to be 
performed with pileup. 
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